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Abstract: For two interacting particles (TIP) in one-dimensional random potential the 
dependence of the Breit-Wigner width T, the local density of states and the TIP localiza- 
tion length on system parameters is determined analytically . The theoretical predictions 
for r are confirmed by numerical simulations. 
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Recently, the problem of two interacting particles (TIP) in a random potential has 
attracted interest of different groups g [||, ||, ||, Q]. It has been shown that two re- 
pulsive/attracting particles can propagate together on a distance l c much larger than 
one-particle localization length l\ in absence of interaction. The first analytical studies 
|l], for TIP with on site interaction on a one-dimensional one channel lattice gave the 
following estimate l c /h ~ Tp ~ (U/V) 2 li, where U is strength of the interaction, V is 
intersite hopping matrix element, p ~ l 2 /V is density of the two-particle states coupled by 
the interaction, and T ~ U 2 /Vl\ is the interaction induced transition rate between these 
states. The numerical investigations [||, |J] definitely confirmed existence of the strong 
enhancement of l c due to interaction. However, a direct verification of the above estimate 
is quite difficult even for the modern computer facilities due to the strong increase of 
required basis with l\. Also the recent numerical results of von Oppen et. al. § and 
Weinmann and Pichard || indicate in one-dimensional case almost linear growth of the 
enhancement factor for l c with U instead of expected U 2 . Due to all these things it would 
be important to have a more rigorous derivation of the factor l c /l\ for this on a first 
glance quite simple problem, at least in a one-dimensional case. To reach this aim we 
started from the computation of the rate V which also characterizes the spread width of 
the Breit-Wigner distribution for eigenfunctions in the basis of eigenstates of noninteract- 
ing particles || |9j. |lOfl . If the parameter dependence of V is known then the ratio l c /l\ can 
be determined from the relation l c /h ~ Tp which have been checked in models of super- 
imposed band random matrices jl], || [9|, |XQ|] . In the present work for calculation of T we 



use the technique developed in 11] which allows to account all orders in the interaction. 
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We consider one dimensional Hubbard model with Hamiltonian 
,t „ , „t „ . \ , ttSTJ J 



na n 

Here a) n is a creation operator of the particle at the site n, V is hopping matrix element, 
and U is on site interaction. We assume that particles are distinguishable and denote the 
type of particle by spin a = ±1/2. Single particle eigenstate is plane wave \p) = ^e ip ™ 
with dispersion e p = — 2V cosp, —it < p < it. We set lattice spacing equal to unity. The 
size of the lattice is denoted by L. 

The Breit-Wigner width can be found in the following way. Forward scattering am- 
plitude / for particles with different spins is given by series of diagrams presented at 
Fig.l. Solid line represents a particle, and wavy line is matrix element of the interaction 
(P3Pa\U\piP2) = x^Pi+P2,P3+P4- Due to optical theorem width of the state \p1P2) = \pi)\P2) 
is related to the forward scattering amplitude: 

r/2 = -Im f. (2) 

One can easily check the coefficient in this relation considering diagram Fig. lb which 
gives usual Fermi golden rule: 

I „ 2 

r— 2 /m/» - -2ME L = (3) 

p 3 p 4 J 4 

„ I - 2 

= 2vr 2^ |(p3P4|LT|p lP2 ) S(E - e 3 - e 4 ). 

P3P4 

Here E is energy of the initial state E = e\ + €2- 

Born term in the amplitude / is given by Fig. la and equals f\ a = U/L. Calculation 
of the diagram Fig. lb is also straightforward 



{P3Pi\U\piP2) 



U 
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f lb X] w.- - M ,n t?Y1 



,. /•' <:i ' 1 • '0 ^ E + 2Vcosp 3 + 2Vcos(p-p 3 ) 



where p = P1+P2 = P3 +Pa is total quasi-momentum. Higher orders in Fig. 1 correspond 
to simple iterations of the box Fig. lb. Therefore summation of the ladder is reduced to 
geometrical progression and the result is 

f{E,P) = l-U/y/&-16V*cos*p/2' (5) 
The scattering amplitude depends only on total energy — 4V < E < AV and total mo- 
mentum — 7r < p < ir. The branch of square root should be chosen in such a way that 
Im f < 0. 

With amplitude (||) one can easily calculate the Breit-Wigner width using optical 
theorem @. But we are interested in the average width at given energy. So we have to 
average over momentum p. Density of the two particle states is of the form 
™ Ldpi [* Ldp 2 



— — / — — 5(p-p 1 -p 2 )5(E + 2V cos pi + 2V cos p 2 ) 

L 2 /{8tt 2 V) 
Vcos^p/2- E 2 /WV 2 ' 

It is nonzero only if square root is real. After integration over momenta we find 



(6) 



The integral in (0) can not be exactly expressed in terms of elementary functions. Pre- 
sented approximate formula is valid with accuracy better than 1% in the interval —4V < 
E < AV . Now we can find the average Breit-Wigner width. 



Y(E) = -2 1m J f(E,p)p{E,p)^-J ' p(E) 

8Vu 2 /L 



(ln4/e + 0.18e) ^\u 2 - e 2 |(l + u 2 - e 2 



F(Z). (8) 



Here u = \U/4V\ and e = |i?/4V| is interaction and energy expressed in units of band 
width: < e < 1. The function F(Z) is defined by 

arctan Z, for u > e 
F(Z) = { ~ (9) 

iln^f, foru<e 



/ | M 2_ e 2| (1 _ e 2) 

V e 2 (l + u 2 -e 2 ) " 
At small energy (e 2 <C u 2 , 1) formula (|8|) can be substantially simplified 
47rT/ 1 u 

r ^~^'mlA' Vl + n 2 ' (10) 
so that at small interaction (e 2 <C u 2 -C 1) it is linear in the interaction. In other 
limit (u 2 <C e 2 , (1 — e 2 )) the width (|8|) is quadratic in the interaction with logarithmic 
correction: 



t, 8V 1 u , 2e\/T^e 2 

T« ; --—In — (11) 

L (ln4/e + 0.18e) e u y ' 

The value of V in ( JTOj ) is significantly larger than in (11) due to the growth of two-particle 
density of states (^) near the center of the band. 

If we now add to the Hamiltonian (|l]) a single particle random potential H ran ^ = 
Wn^n(T a na with a disorder homogeneously distributed in the interval — W < w n < W , 
then one particle eigenstates in infinite lattice become localized with localization length 



l\ ~ 24(V/W) 2 yl — ef/AV 2 , where e\ is one particle energy. However as soon as l\ S> 1 
the above calculation of the average width remains valid. The reason for this is that 
l\ S> 1 is the only condition which we need to formulate scattering problem and to use 
conventional diagram technique. Distribution of V depends on the relation between size 
of the box L and the localization length l\. If L < 1% all values of T are of the order 



of the average value given by (|8|). For L 3> h the average value is still given by (||). 
However in this case T vanishes for majority of the states . These are the states in which 
particles are localized far from each other and practically do not interact. On other hand 
the width for the states with interparticle distance of the order l\ is approximately the 
same as for particles in a box of size L ~ l\ so that V is given by eqs.(||),(|9|) with L 
replaced by l\. The two-particle localization length l c for such states is determined by 
the relation l c /h ~ T(E)p(E), with T calculated at L ~ l\. This relation is valid if 
many unperturbed states are mixed by interaction |2| so that T{E)p(E) > 1. In the 
opposite case T(E)p(E) <C 1 the above relation is not valid (8|, ^, 10 1 and the interaction 



can be treated in a perturbative way. In this regime "Rabi oscillations" between two 
quasi-degenerate levels play an important role ||. 

Above we have considered distinguishable particles. The generalization to identical 
particles is rather simple: the width T vanishes if coordinate wave function is antisym- 
metric, and it is doubled in comparison with eqs.®, @ if coordinate wave function is 
symmetric. 

To check the above theoretical formula for the Breit-Wigner width T we studied 
numerically the model (|l|) of two identical interacting particles (symmetric coordinate 
wave function) in the disordered potential on a ring of size L which is less or comparable 
with one-particle localization length l\ ~ 24(V/W) 2 . Using Lanczos technique (see for 
example [^]) we determined the local density of states for symmetric configurations in 
the basis of noninteracting eigenstates: 



p w (E-e mi -e m2 ) = J2\Mmum 2 )\ 2 S(E - E x ) (12) 

A 



Here E\ is the eigenenergy of TIP while e mi 2 are one-particle eigenenergies. The depen- 
dence of pw on E is well described by the Breit-Wigner distribution 

Pw{E) = 2^ + r74] (13) 

an example of which is shown in Fig. 2 . The comparison of numerically obtained T with 
theoretical prediction @, (||) in the regime T(E)p(E) > 1 is shown in Figs. 3, 4 for different 
energies as the function of interaction. The theory gives good agreement with numerical 
results for 15 < L < 300 and variation of scaled width TL/V by more than 2 orders 
of magnitude. For the states with the energy close to the band center (E ~ 0) (Fig. 3) 
the dependence of Tp on U is almost linear for U < V (see (Q), (|l~0|)). Therefore, the 
TIP localization length l c according to the relation l c /h = CYp ~ 2Ch(U/V)/ir also 
varies linearly with U. Here, we took the values of T and p at L = 1% and introduced the 
numerical coefficient C to take into account the uncertainty of this choice. According to 
the numerical results at the center of the band l c /h ~ 0.2li(U/V) which is in good 
agreement with the above theoretical expression and gives C ~ 1/4. 

For energies away from the band center and small interaction \U\ <ti \E\ the enhance- 
ment factor according to 0, ([□]) is l c /h w hU 2 \n(2E /U) / (Att 2 V E) where we have used 
the above value of C. The dependence on U is almost quadratic in agreement with the 
first estimate [jl|, Q. However, due to the logarithmic correction, to observe clearly the U 2 
behavior one should go to really small U values and since the condition Tp > 1 should be 
also satisfied this can be reached only for quite large values of l\ or L. In this respect our 
numerical approach based on the measurement of T is more efficient than the one used in 
Q . It allows to see the behavior U 2 In U away from the band center in agreement with 
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the theory (||), (^|) (see insert in Fig. 4). At moderate U/V > 0.3 values in the presence of 
numerical fluctuations the dependence of T on U is hardly distinguishable from a linear 
one (see normal scale in Fig. 4). In our opinion this is the reason why the linear behavior 
in U had been attributed in [|j also to the states away from the band center. As for 
the result of Ref. Q the system size there was too small (L = 25) and the main part of 
the data (Fig. 4 with U/V < 0.4) corresponds to the different regime Tp < 1. In this 
perturbative case the typical energy scale which determines the change in level statistics 
is determined by Rabi oscillation frequency in a pair of quasi-degenerate states which is 
proportional to U ^ . Also, one should keep in mind that the results there are integrated 
over the whole energy band including the center of the band where the dependence on U 
is linear even for Tp > 1. 

Turning back to our numerical data (Fig. 4) we would like to mention that there is a 
significant difference from the theory for negative U < —1.. Generally, we should expect 
such difference for \U\ S> 1 when the spectrum is composed from two separated energy 
bands and the basis of plane waves used for computation of width T becomes inadequate. 
For example, in this regime the density of states is not described by (^). However, we 
cannot say why this change goes in so asymmetric way for negative and positive U while 
for \U\ < 1 the width T is independent of sign U in agreement with the theory. We would 
like to note that such asymmetry for attraction and repulsion away from the band center 
and relatively strong interaction U ~ V has been seen recently in [|| for the ratio l c /h- 
Also a change in the behavior of T has been observed in || for U > V. 

In summary, taking diagrammatically into account the effects of interaction we have 
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derived the analytical formula for the Breit-Wigner width T which determines the en- 
hancement factor l c /h ~ > 1 for TIP in one-dimensional random potential. Our 
analytical and numerical approaches can be also used for calculation of the TIP width in 
2- and 3-dimensional disordered systems where according to Imry estimate [§] interaction 
between two quasi-particles can strongly affect transport properties. 
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Figure captions 

Fig. 1: Diagrams for the forward scattering amplitude / in (2) - (5). 

Fig. 2: Local spectral density pw(E) computed for the TIP eigenstates in the energy 
interval [-0.1, 0.1] for the case L = 150, U = 1, V = 1, and W = 0.4. The full line 
gives the best Breit-Wigner fit (13) with T = 0.0073. The theoretical prediction is 
T = 0.0072. 

Fig. 3: Scaled Breit-Wigner width TL/V as a function of the rescaled interaction 

computed in the energy interval E/V € [—0.1,0.1]. The system size is L = 15 
(W/V = 1., empty circles), L = 25 (W/V = 1., empty squares), L = 40 (W/V = 
0.6, empty diamonds), L = 60 (W/V = 0.5, full circles), L = 80 (W/V = 0.5, full 
squares), L = 100 (W/V = 0.5, full diamonds) L = 150 (W/V = 0.4, full triangles 
up) and L = 200 (W/V = 0.35, full triangles down). The solid line gives the 
theoretical prediction (8), (9) multiplied by 2 to take symmetrization into account. 

Fig. 4: Scaled Breit-Wigner width TL/V as a function of the rescaled absolute value of 
the interaction ^ computed in the energy interval E/V G [1., 1.2]. The system 
size is L = 15 (W/V = 1., empty circles), L = 25 (W/V = 1., empty squares and 
W/V = 0.5, empty diamonds), L = 40 (W/V = 0.8, empty triangles up and W/V = 
0.5, empty triangles down), L = 60 (W/V = 0.5, full circles), L = 80 (W/V = 0.5, 
full squares), L = 100 (W/V = 0.5, full diamonds), L = 100 (W/V = 0.5, negative 
U, crosses), L = 150 (W/V = 0.4, full triangles up), L = 200 (W/V = 0.35, full 
triangles down) and L = 300 (W/V = 0.25, full triangles left). 
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